Fringe free holographic measurements of large amplitude vibrations 
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In the measurement of the ampUtude of vibration of objects, holographic imaging techniques 
usually involve fringe counting; due to the limited resolution of images, measurements of large 
amplitudes are not accessible. We demonstrate a technique that suppresses the necessity of fringe 
counting: frequency sideband imaging, where the order of the sideband is considered as a marker of 
the amplitude. The measurement is completely local: no comparison with another reference point 
on the object is necessary. It involves a sharp variation of a signal, which makes it robust against 
perturbations. The method is demonstrated in an experiment made with a vibrating clarinet reed; 
phase modulations as large as 1000 radians have been measured. 
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The observation of interference fringes in holographic 
methods provides accurate measurements of the ampli- 
tude of vibration of objects. Powell and Stetson [1] have 
shown that the fringes on the holographic reconstruction 
of a vibrating object correspond, after time averaging, to 
zeros of the Bessel function Jo(^), where (I>{x,y) is the 
amplitude of phase modulation of the optical field emit- 
ted by the object at point x,y. Digital holography was 
introduced in 1994 by Schnars and Jiipner [2]. In 2003 
Picard et al. [3] transposed the time averaging method 
to digital holography. 

In a previous letter [4], we described sideband digi- 
tal holography, based on the detection of the light back- 
scattered by a vibrating object at different sideband fre- 
quencies; the fringes for sideband n then correspond to 
the zeros of the n-th order Bessel function J„(<?). As 
in the work made by Aleksoff [5] in 1971, the reference 
beam was frequency shifted to select one sideband n, but 
the use of acousto-optic modulators and numerical tech- 
niques provided much more flexibility. In reference [4], 
we have shown how the comparison of dark fringes for dif- 
ferent sideband leads to a determination of the vibration 
amplitude <l>{x,y) at each point of the object. This de- 
termination is non-local, since it involves counting fringes 
from one reference point of the image to the point of in- 
terest, so that large amplitudes are not accessible. 

In this letter, we demonstrate another approach that 
completely eliminates the necessity of counting fringes; it 
gives a local measurement of the amplitude of vibration, 
even for large values. For any pixel of coordinate x, y, we 
consider the sideband order n as a variable and we plot 
the intensity / as a function of n. One can then easily 
determine <P since I{n) exhibits a sharp variation from 
maximum to zero near n ~ <!>. The method is robust and 
can easily be used even when the fringes become so nar- 
row that they cannot be resolved, which gives immediate 
access to large amplitudes of vibration. 

Consider a point of the object vibrating at frequency 



VA and amplitude z^^y,; its displacement z(t) is: 

= z„,,,^ sin(27ri//it) (1) 

In backscattering geometry, this corresponds to a phase 
modulation: 

(p{t) = 47rz(i)/A = ^sin(27ri^At) (2) 

where A is the optical wavelength and ^ = Anz^-,^,-/ \. 
The scattered field is then: 



E{t) = £:e^>»*+'^(*)) = f ^ J„(<?)e 



j[vo+nvA)t 



(3) 



where £ is the complex amplitude of the field, the 
frequency of the illumination optical field, and J„ the 
rt-th order Bessel function of the first kind; J_„(z) = 
(—1)" Jn{z) for integer n and real z. The scattered field 
is then the sum of sidebands with frequencies z/q + nuA 
and intensities /„ given by: 



/„(<!>) = \8JnmY 



(4) 



Figure 1 (a) shows the intensity of the sidebands as a 
function of n, assuming <!> ~ 30.3 rad. When n is consid- 
ered as a continuous variable ii = [v — vq)Iva giving the 
Doppler frequency shift in units of VA^ one obtains the 
light grey shade corresponding to the Doppler spectrum 
deduced from the vibration velocity distribution. This 
continuous spectrum is proportional to [ 1 — (n/^)^]~^/^, 
where n is confined between the values that corre- 
spond to the maximum Doppler shift. The discrete spec- 
trum has a similar behavior, remaining mostly confined 
between the same n values, and dropping abruptly from 
a maximum reached close to n = ±<P to almost zero. 
This is the key idea of our method: we measure the fre- 
quency position of this sharp variation and deduce from 
it the value of the vibration amplitude. This method can 
be seen as a discrete spectrum version of laser Doppler 
imaging of non-periodic motions reported in [6, 7]. 
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FIG. 1: (a) Relative intensities of the sidebands as a function 
of n for fixed $ — 30.3 rad. The vertical lines show the 
intensities of the discrete n components of the real spectrum. 
The light grey shade shows the Doppler spectrum obtained 
from the vibration velocity distribution, with a continuous 
variable on the horizontal axis n = — i'o)/i^A- Both spectra 
fall abruptly beyond n — ±30.3, which corresponds to the 
Doppler shift associated with the maximum velocity, (b) The 
broken line shows the values ni as a function of <1>, where ni 
is the value of n giving the maximum intensity in the discrete 
spectrum; the full line shows n2, where 722 is the value of n 
for which the intensity if half the maximum, which gives a 
very good approximation of (c) The broken and full lines 
respectively show ni — $ and n2 — ^ as a function of 



Figs. l(b,c) give more detail on this sharp variation. 
For each value of we calculate the value ni(^) of n 
that, in Eq. (4), gives the maximum intensity. To avoid 
the quantization noise induced by discrete variables, the 
calculation is made with Bessel functions of fractional or- 
der, but of course only integer values of n are relevant to 
the experiment. It is known [8] that ni(^) ~ within 
a correction of order ; the broken lines shows ni(<?) 

(b) , and the difference iii{<P) — <P (c). To determine more 
precisely the location of the abrupt drop from maximum 
to zero, we calculate the n value n2(^) for which the 
intensity is half the maximum. The 712 {'P) and the differ- 
ence 77,2 (^) — ^ are shown by the full line in Fig. 1 (b) and 

(c) respectively. Clearly, 7x2 gives an accurate evaluation 
of ^. 

The experimental setup is the same as described in 
Ref. [4]. As in [4], we have chosen a clarinet reed as 
the vibrating object to experimentally demonstrate the 
method. The reed, vibrating on its first flexural reso- 
nance mode at frequency va ~ 2100 Hz, is illuminated at 
A = 650 nm by a laser field Ej at frequency i/q. The CCD 
camera (frame frequency vccd = 10 Hz) records the in- 
terference pattern (the hologram) between the scattered 
light and the local oscillator beam (field i?Lo, frequency 
^Lo)- Two acousto-optic modulators (Bragg cells) with 
frequencies vaomi and vaom2 are used to adjust the fre- 
quencies vq ^ VL + VAOM2 and vlo = vl + vaomi] an 
arbitrary sideband n can then be selected by adjusting 
these frequencies. A 4-phase detection of sideband n is 
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FIG. 2; Cube of data obtained from the reconstructed holo- 
graphic images of a vibrating clarinet reed; sideband images 
with n = 0, 20, 40. ..120 are shown in arbitrary linear scale. 
By choosing n, one moves the border of the illuminated re- 
gion on the object, obtaining a local marker of the amplitude 
of vibration. The white dashed lines correspond to x = 249 
and y = 750, i.e. to the point chosen for Fig. 3. 

obtained by adjusting them to fulfil the relation: 

+ nvA - VLoin) = vccdI^ (5) 

We record a sequence of consecutive CCD images 
lo, Ii, ■■■Im-1, where M is a multiple of 4. The complex 
hologram H in the plane of the CCD is is then: 

M-l 

H=j2r im (6) 

m=0 

From H, the images are reconstructed as in Ref. [4] by 
the standard convolution method [2] , which provides the 
map of the complex field E{n, x, y) in the plane of the ob- 
ject. By successively adjusting the frequency VLoin) of 
the local oscillator to appropriate values, we then record 
the intensity images \E{x,y,n)\^ of the sidebands as a 
function of x, y and n. We then obtain a cube of data 
with three axes x, y and n, where x, y are expressed in 
units of pixels of the reconstructed image of the reed. Fig- 
ure 2 shows the images obtained for for n = 0, 20, ...120 
that correspond to cuts of the cube along x,y planes. 
The right part of the reed {x > 800) is clamped on 
the mouthpiece. The images illustrate how, when n in- 
creases, the fringes move towards regions with larger am- 
plitudes of vibrations: no signal is obtained in regions 
where = Attz^^^/X < n. This well known property of 
Bessel functions allows one to get a marker on the ob- 
ject, signalling regions where the amplitude correspond 
to <P ^ 120 rad. Fig. 3(a) displays a cut of the cube of 
data along the x, n horizontal plane y = 750 (horizontal 
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FIG. 3: Images corresponding to cuts of 3D data of the recon- 
structed images along the planes y = 750 (a) and x = 249 (b). 

Fig (a) shows the deformation of the object along its axis and 
Fig (b) a transverse cut with a slight vibration asymmetry. A 
logarithmic intensity scale is used. 




FIG. 4: (a): image reconstructed with sideband n = 330, 
with a large amplitude of vibration, (b) is the equivalent 
of Fig. 3(a), but with positive and negative n values; one 
measures a maximum vibration amplitude of z„ax — 60 yum. 
A logarithmic intensity scale is used. 

white dashed hnc in Fig. 2). The envelope of the non-zero 
(non black) part of the image provides a measurement of 
the amplitude of vibration in units of A/47r. We actually 
obtain a direct visualization of the shape of the reed at 
maximal elongation, from the right part clamped on the 
mouthpiece to the tip on the left. The maximum am- 
plitude correspond to <P ~ 120 rad. Fig. 3(b) displays a 
transverse cut along the y, n vertical x = 249 plane (ver- 
tical white dashed line in Fig. 2); a slight asymmetry of 
the reed vibration is clearly visible. 

Figures 2 and 3 have been obtained by exciting the 



reed at the resonant frequency ua = 2123 Hz and at 
moderate excitation level. The consistency of amplitude 
measurements obtained in this way with those obtained 
by fringe counting is illustrated in Fig. 3 of Rcf. [4]. 
We have also used higher excitation amplitudes, about 
10 times larger (the frequency of reed resonance is then 
shifted to i^a = 2020 Hz). Figure 4 (a) shows the images 
obtained for n = 330: the fringes are now completely 
unresolved, but the transition from zero to non-zero in- 
tensity remains very clear. With a single hologram, and 
without fringe counting, one obtains a clear marker of 
the line where ^{x, y) = 330 rad. 

Fig. 4(b) shows the equivalent of Fig. 3(a), but with 
a higher excitation level, and this time for positive and 
negative values of n; data range up to about \n\ ~ 1140, 
corresponding to z^^^ ~ 58.4 /im. Since the vibration 
amplitude is much larger than A, the continuous approx- 
imation for n is valid, and the images of Fig. 3 can 
be reinterpreted in term of classical Doppler effect: if the 
signal intensity In{x, y) at frequency vq+tipa is non-zero, 
at some time t of the periodic motion the Doppler shift 

+ 2v{x, y,t)/c is equal to Wa {v and c are the reed and 
light velocities). The accuracy of the measurement of the 
amplitude corresponds to about 0.5^ (see discussion in 
the supplementary material). A fine analysis of the 3D 
data shows a slight asymmetry between the positive and 
negative n values. In Fig. 4(b) for example, the maxi- 
mum value for \n\ is about 1160 for n > 0, and 1120 for 
n < 0. This means that the motion of the reed in not 
perfectly sinusoidal, and that maximum velocities (i.e. 
maximum Doppler shifts) are slightly different for the up 
and down motion of the reed. A more detailed study of 
the modification of the spectrum due to the superposi- 
tion of several harmonic motions can be done, but falls 
beyond the scope of this letter. 

In conclusion, taking advantage of the frequency (or 
the sideband order n) of the light scattered by a vibrat- 
ing object adds a new dimension to digital holography. 
Each pixel of the image then provides an information 
that is completely independent of the others, which re- 
sults in redundancy and robustness of the measurements. 
Looking at the edges of the spectrum provides an accu- 
rate determination of the vibration amplitude and avoids 
a cumbersome analysis of the whole cube of data cube, 
giving easy access to a measurements of large amplitudes 
of oscillation. Our method could be combined with other 
techniques, such as phase modulation of the reference 
beam [9, 10], in order to get information on the phase of 
the mechanical vibration without loosing the robustness 
of the measurement. 
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